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We report a method to obtain the stress of crystaUine materials directly from lattice deformation by Hooke's 
law. The lattice deformation was calculated using the crystallographic orientations obtained from electron 
backscatter diffraction (EBSD) technology. The stress distribution over a large area was obtained efficiently 
and accurately using this method. Wurtzite structure gallium nitride (GaN) crystal was used as the example 
of a hexagonal crystal system. With this method, the stress distribution of a GaN crystal was obtained. 
Raman spectroscopy was used to verifj^ the stress distribution. The cause of the stress distribution found in 
the GaN crystal was discussed from theoretical analysis and EBSD data. Other properties related to lattice 
deformation, such as piezoelectricity, can also be analyzed by this novel approach based on EBSD data. 

Stress intensity and distribution of crystalline materials seriously affect their properties and the performance 
of devices based upon them. For example, the bandgap of a compound semiconductor crystal, which 
determines the emission wavelength of devices, is affected by the stress \ Different stress conditions, 
compressive or tensile, result in different properties. Lattice deformation, represented by changes in the lattice 
parameters, leads to the introduction of stress^'^. The stress condition is also important for adjusting crystal 
growth conditions at different growth stages. So it is important to obtain the stress information of crystalline 
materials. Various methods have been used to obtain the stress information of crystalline materials. The lattice 
parameters of crystals can be detected by the X-ray diffraction (XRD) method directly from the angles of 
diffraction peaks^. The stress can be calculated from changes in the lattice parameters between the stressed 
and stress-free lattice. However, because the radius of the incident X-ray is just several microns, only a small 
area can be tested by this method. As such, it is difficult to obtain a large area mapping of stress in crystalline 
materials by XRD. The X-ray microdiffraction method increases the spatial resolution and can be employed to 
analyze the stress of crystalline materials^. High resolution transmission electron microscopy (HR-TEM) can also 
be employed to analyze lattice parameters directly^. However, Just like XRD, it is also difficult to obtain stress 
distribution information over a large area using HR-TEM. It is also especially challenging to prepare crystalline 
samples for analysis by HR-TEM. At the same time, several indirect methods have been exploited to analyze the 
stress condition of crystalline materials. For example, Raman spectroscopy has been used to evaluate the stress of 
crystalline III-V materials, like gallium nitride (GaN), utilizing the Raman shift wave number distinction of 
E2(high) mode between stress and stress-free GaN crystals^. However, the stress value obtained depends on the 
measurement conditions and the selection of the stress-free GaN crystal E2(high) peak position. It is difficult to 
identify a stress-free crystal accurately. Photoluminescence^ and micro -reflectance spectroscopy^ have also been 
employed to analyze the stress of crystalline materials. Stress values obtained by these two methods are generally 
similar to those obtained from Raman spectroscopy. However, the stress values obtained from these three indirect 
methods are not acquired from the lattice parameters directly. 

The electron backscatter diffraction (EBSD) technique is an effective method to obtain information concerning 
crystal orientation and for phase identification^°'^\ EBSD is also very useful for the microstructural character- 
ization of crystalline materials. Because of the high spatial resolution and good strain sensitivity, EBSD has been 
used to measure the crystal orientation in compound semiconductor hetero- epitaxial structures^^'^^. The distri- 
bution of the crystal orientation over a large area can be obtained easily and quickly from EBSD mapping. The 
orientation data at each point serves as the foundation for analysis of the stress distribution in crystalline 
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materials. With this method, the stress value is obtained directly from 
the crystal structure because the orientation data given by the geo- 
metry of EBSD Kikuchi patterns depends on the crystal structure. 

In this article, we report a method for analysis the stress distri- 
bution of crystalline materials directly from the lattice deformation 
derived from the crystallographic orientation identified by the EBSD 
technique. The method offers a potential to analyze the properties 
associated with lattice strain. Hexagonal GaN crystal grown by 
hydride vapor phase epitaxy (HVPE) on a foreign substrate was used 
as the study sample. The stress distribution of GaN crystal was 
determined according to the degree of misorientation at different 
positions of the hetero- epitaxial GaN crystal. The stress conditions 
obtained by Raman spectroscopy verified the results calculated from 
the EBSD measurements. 

Experimental 

An EBSD system (Oxford Instruments INCA Crystal EBSD system, 
Nordlys EBSD Detector and HKL CHANNELS software) mounted 
on a field emission scanning electron microscope (FE-SEM, Hitachi 
S-4800) was used to obtain the strain information of HVPE-grown 
GaN. The EBSD measurement was operated at 20 kV with a work 
distance of 20 mm and a sample tilt of 70°. 

MOCVD grown GaN on a (0001) sapphire wafer was employed as 
the substrate for the HVPE process. Patterns on the surface of the 
MOCVD-GaN were hexagonal mask arrays with diameters of 
5.4 |im and a center-to-center distance of 8.1 jim. The patterned 
substrate was set in a vertical HVPE system with bottom-fed gas after 
a cleaning process. The thickness of the GaN crystal obtained by 
HVPE was 430 |im. After a cooling process, the GaN crystal self- 
separated from the substrate. 

Cross-sectional mapping of the free standing HVPE grown in the 
<0001> direction was conducted by EBSD. The mapping area was 
60 |im X 260 |im employing a step length of 1 |im. Z axis mapping 
by Raman spectroscopy was also employed to identify the stress 
distribution of the GaN crystal. 

Results 

EBSD mapping of a cross-section of GaN crystal was obtained, as 
shown in Fig. 1 (a). The <0001> direction was set as the reference 
direction. Deviation from the <0001> direction was obtained from 
the spatial orientation information collected from backscattered 
Kikuchi diffraction patterns (three Euler angles) at each EBSD map- 
ping point. The EBSD results show that the deviation from the 
<0001> direction near the interface is larger than that far from it. 
The mapping color near the interface is yellow and red, which indi- 
cates a deviated from the <0001> direction of 6.2°. As the distance 
from the interface increases, the mapping color changes from yellow 
to green and the deviation from the <0001> direction reduces to 
5.8°. This result demonstrates that the extent of misorientation near 
the interface is greater than that far from the interface and deviation 
from the <0001> direction decreases with increasing distance from 
the interface. 

Discussion 

In crystalline materials, mechanical stress and deformation are, 
within the limit of Hooke's law, proportional to one another. For 
sufficiently small stress and deformation, the relationship of these 
quantities can be described, when infinitesimal torques or rotations 
are neglected^^, by 

^ij = ^ijki^ki (1) 

In equation (1), c^ki is the elasticity tensor whereby the stress tensor 
(Gij) can be calculated from the deformation tensor (ski). The 
deformation is described by changes in the lattice parameters 
between the stressed and stress-free lattice. The elasticity tensor is 
a fourth-rank tensor with 3^ = 81 components, although symmetry 
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Figure 1 | Crystallographic orientations represented as deviations from 
the <0001> direction obtained from EBSD mapping data (a), and the 
calculated stress values based upon the EBSD crystallographic orientation 
mapping results (b). 

reduction decreases the number of independent components.The 
number of independent components and, hence, the expression of 
Cijki vary according to different crystallographic symmetry condi- 
tions. The stress can be calculated by this relationship using a mater- 
ial-specific elasticity tensor. 

Single crystal GaN, which is a wurtzite structure crystal, possesses 
Cev^ space-group symmetry, and is an example of a hexagonal crystal 
system. The stress distribution of GaN crystal on a foreign substrate 
is apparent. The GaN crystal was grown on a sapphire substrate along 
the <0001> direction. However, the actual growth direction 
deviated slightly from the <0001> direction due to the mismatch 
of the lattices and thermal expansion coefficients between GaN and 
the sapphire substrate. Stress was introduced because of the resulting 
misorientation^^'^^. Dislocations and cracks released some of the 
stress local to the defect, but residual stress remained in the GaN 
crystal. The distribution of the stress along the <0001> direction 
indicate that the stress decreased with increasing distance from the 
original interface between the GaN crystal and the substrate (N polar 
face), and the smallest value of stress was observed at the top free 
surface of the GaN (Ga polar face)^. The stress distribution is distinct 
in this kind of hetero -epitaxial crystal. In this case, EBSD mapping is 
appropriate for stress calculation of the GaN crystal grown on a 
foreign substrate. 

EBSD mapping was used to determine the crystallographic ori- 
entation at regular points on the sample for calculation of the stress. 
In EBSD mapping, the crystallographic orientation is recorded using 
three Euler angles that describe a minimum set of rotations that can 
bring one orientation into coincidence with another. During a mea- 
surement, this is the relationship between the EBSD detector and the 
particular point on the sample being measured^^. Bunge convention 
is most commonly used to describe the Euler angle rotations. The 
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three Euler angles ((pi, O and 92) represent the following rotations, 
which are shown schematically in Fig. 2(a): first, a rotation of (pi 
around the Z-axis; second, a rotation of O around the rotated X-axis; 
and third, a rotation of 92 around the rotated Z-axis. The ideal Euler 
angle values for a cross -section of a GaN crystal grown on a sapphire 
substrate along the <0001 > direction are (pi = 180°, O = 90° and (p2 
= 0° as shown in Fig. 2(b). In this case the cross-section surface of the 
GaN crystal is identified as (11-20) in the hexagonal crystal system. 
However, the actual crystallographic orientation of GaN deviates 
from the ideal value, as shown in Fig. 2(c). 

Misorientation of the GaN crystal was determined by the differ- 
ences between the ideal Euler angles and the actual values. The EBSD 
mapping data show that the three Euler angles of GaN deviated from 
the ideal values. These deviations are given as: 

A(pi = 180°-(pi, AO = 90°-<I), A(p2 = (p2. 



The angle deviations indicate the degree of misorientation, and 
they affected the lattice parameter at the projections on the (0001) 
plane (Fig. 2 (c)). The lattice parameters of the projection can be 
calculated from the rotation matrix. The relationship between the 
Bunge-Euler angles and the rotation matrix^^ is 



cos ^1 COS (P2 — sin cos ^ sin (p 
sin (p^ cos (P2 + cos cos ^ sin (p 
sin ^ sin ^2 



- cos (p^ sin (P2 — sin cp^ cos ^ cos ^2 sin sin 

- sin sin ^2 + '^os cos ^ cos ^2 — sin^cos^^ I ' (2) 



which can be described in terms of the misorientation angles as: 

- cos A(9i COS ^(p2 — sin A^j sin Ai^ sin A^2 '^^^ ^Vi sin A^j ^ sin A.(p^ sin A(/> cos A^2 '^os A(/> sin A^j \ 
^ = I sin A^?! cos ^(p2 — cos A^j sin A(/i sin A^j ^ sin A^^j sin A — cos Acp^ sin A(/i cos A^2 cos Ac/i cos /^cp^ j • ( 3 ) 
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Figure 2 | Euler angle rotations according to Bunge's convention (a), and 
ideal Euler angle values of a cross-section of free-standing GaN crystal 
self- separated from the sapphire substrate along the <0001> direction 
(b). Comparison of the lattice projections on the (0001) plane for the ideal 
and actual crystallographic orientations (c). 



The lattice parameter is determined as the distance from the origin 
point to point P ( — ^ ao, y ,0) in the initial coordinates. Based upon 
the ideal condition ((pi = 180°, 0) = 90° and (p2 = 0°), the P rotates to 
point P' { — y^ao,0,y), indicating that the c-plane moves to the 
XOZ plane in the initial coordinates. Taking the misorientation into 
account (using equation (3)) point P rotates to point P" according to 
the following: 

P'=g-P\ (4) 
where T represents the transpose, resulting in 

/ |cos A9Jj( sin +'v/3 cos A^jj)— f sin sin cos ~ "v/^ sin A^2) \ 

P"^ = — ^^(sinA^j cos A^2 ~ <^os A^j sinA(^sinA9J2)— f (sinA^j sinA^2 + '^os A^jj sinA^cos A^2) ' (5) 
y — ^ cos A^ sin A^2 + f ^'^ '-"^ ^Vi / 

The XOZ plane is set as the basis plane, and in this case, using the 
projection of P" on the XOZ plane, the lattice parameter is 

a=—{ [cos A^i ( sin A^2 + ^^2) ~ sin A^i sin A(^( cos A^2 



-V^sin A^2)l^ + [cos A(^(a/3 sin A^2 ~ cos A^2)l^}^ 



(6) 



where a is the actual lattice parameter and ag is ideal value. 

To the space-group symmetry of single crystal GaN, the strain 
tensor e is diagonal and possesses the components: 



flo)/flo, 
(c-Co)/co. 



(7) 
(8) 



In these equations, a and c are the lattice parameters of single 
crystal GaN, where ag and Cq are the lattice parameters in a condition 
of equilibrium. The corresponding diagonal stress tensor compo- 
nents are given by Hooke's law in the limit of small deviations from 
equilibrium as 



— <Jyy — (Cii + Ci2)e^^ + ClsS^; 



(9) 
(10) 



Four of the five independent stiffness constants Qj of the consid- 
ered wurtzite crystal using in equation (9) and equation (10) are 
shown in Table 1^^'^°. The biaxial stress in the plane perpendicular 
to the c axis of the lattice is described by constant forces in this plane 
such as 

(^xx = (^yy^ (11) 

a,, = 0. (12) 



Then, Hooke's law (equation (9) and equation (10)) provides a 
relationship between the strain components is 



_2Ci3 

^33 



(13) 



Table 1 


Elastic constants (GPa) of wurtzite GaN crystal 


GaN(Wz) 


Cn 


390^9 


C,2 


145^9 


C,3 


10620 


C33 


398^9 
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The biaxial modulus is given in terms of the elastic stiffness constants 
as 



- (Ci 1 + Ci2 — )£x. 

^33 



(14) 



The relationship between the stress and strain is obtained from 
equation (14). 

The stress obtained via equation (7) by EBSD (from equation (6)) 

is 

(Txx = (Ci 1 + Ci2 -^)(a- ao) 

^33 

2(^2 ^ 

= (Ci 1 + Ci2 - ){ 1 ~ 2 ^ ™ ^ ^^^^ (1 ^) 

— sin A(p^ sin A^( cos A^2 ~ ^^2))^ 

+ ( cos A(/>(\/3 sin A^2 ~ cos A^2))^l^}- 

The stress obtained from equation (15) is shown in Fig. 1 (b), the stress 
value is large near the interface, decreases towards the middle region, 
and, approaching the Ga-polar face, the stress becomes small and 
constant. 

Raman spectroscopy was also carried out to evaluate the stress 
distribution in the hetero- epitaxial GaN crystal. For wurtzite struc- 
ture GaN, E2(high) mode of Raman spectroscopy is considered to be 
affected only by stress. A biaxial stress of one GPa shifts the E2(high) 
Raman mode by 4.2 cm~\ and the wave number of the stress-free 
GaN E2(high) mode^^ was taken as 566.2 cm"\ The free-standing 
GaN crystal self- separated from the substrate was measured by 



Raman spectroscopy at different depths from the N-polar face 
towards the Ga-polar face along the <0001> direction. The peak 
position of E2(high) Raman mode at different depths and the cor- 
responding stress are shown in Fig. 3 (b), where a depth of 0 |im 
represents the surface of the N-polar face, and increasing distance 
indicates motion towards the Ga-polar face. The peak position of 
E2(high) Raman mode near the interface (0-70 jim) exhibits a higher 
wave number than it far from the interface (110-200 jim). These 
results demonstrate that the stress was larger near the interface than 
near the free surface, and a middle transition region (70-110 jam) is 
observed. 

The stress distribution of the GaN crystal obtained from Raman 
spectroscopy is consistent with the results calculated from EBSD 
data. Both the results obtained from these two methods demonstrate 
that the stress near interface is large, obviously decreases over a 
middle range, and maintains a small value over some distance from 
the Ga-polar face. Clearly, the impact of the sapphire substrate on the 
stress of the GaN crystal occurs over a limited range of distance. As 
such, the impact becomes a small and constant for a sufficiently thick 
GaN crystal because the lattice parameter approaches a value close to 
that of the stress-free condition sufficiently far from the interface. 
Therefore, the lattice mismatch described by Equation (6) also 
becomes small. It should be noted that the stress value calculated 
from EBSD is larger than that obtained from Rama. This is because 
the stress value obtained from Raman is a relative value that depends 
on the selection of the E2(high) mode wave number associated with a 
stress free crystal. However, it is difficult to identify a point on the 
crystal surface that is completely stress-free. Therefore, the stress 
determined by empirically selecting the E2(high) mode wave number 
of a stress-free crystal is not absolutely accurate. However, the influ- 
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Figure 3 | The stress values calculated the crystallographic orientations at different positions in the cross-section of the free-standing GaN crystal with 
respect to the distance from the N-polar face (a), and the E2(high) Raman mode peak position at different distances from N polar face and the 
corresponding stress (b). 
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ence of different kinds of stress on the EBSD data is not the same. The 
calculated stress will be underestimated when based on misorienta- 
tion in some condition. Meanwhile, other properties related to lattice 
deformation, such as piezoelectricity, can also be analyzed by this 
approach based on EBSD data. The only difference is that the elasti- 
city tensor will be replaced by some other tensor. The good sensitivity 
makes analysis on a micron scale convenient, at the same time; 
properties distributed over a large area can also be obtained by this 
method. 

Conclusions 

A method for calculating stress in crystalline materials directly using 
the lattice deformation identified by EBSD was provided. The stress 
of crystalline materials at each mapping point was obtained from 
EBSD by this method. The stress distribution over a large area was 
obtained efficiently and accurately. Crystalline GaN of a wurtzite 
structure grown by HVPE on a foreign substrate was used as an 
example of a hexagonal crystal system. The stress distribution of 
the GaN crystal obtained from Raman spectroscopy confirmed the 
results provided by EBSD. We believe that the stress distributions of 
other crystalline materials can also be calculated using this method 
simply by changing the form of the elasticity tensor. 
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